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Abstract. We introduce a class of partial differential equations on metric graphs associated with mixed 
evolution: on some edges we consider diffusion processes, on other ones transport phenomena. This yields 
a system of equations with possibly nonlocal couplings at the boundary. We provide sufficient conditions 
for these to be governed by a contraction semigroup on a Hilbert space naturally associated with the 
system. Subsequently, we devote our interest to further natural questions in the fields of parabolic and 
elliptic differential equations and eventually obtain a characterization of the sub-Markovian property and 
a description of the spectrum. We conclude the article showing that our setting is also adequate to discuss 
specific systems of diffusion equations with boundary delays. 



1. Introduction 

In the literature, usually considered problems concern networks whose ongoing dynamical processes are 
homogeneous: on each link the evolution is modeled as a transport, a diffusion, a wave, a beam, ecc. 
However, many physical models consist of coexisting, interacting processes of different type. On different 
edges a different kind of dynamics may take place; or else, one may introduce fictitious, auxiliary edges in 
the model in order to describe certain phenomena (like delays) in a more efficient way. Accordingly, our 
aim is to discuss a Cauchy problem for a system of partial differential equations of mixed type associated 
with an operator A: a part of the system (denoted by E<j) satisfies a one-dimensional heat equation whereas 
the remaining part (denoted by E t ) satisfies a one-dimensional transport equation, hence 



We assume the interactions of the different subsystems to take place only at the boundary. In this way 
we can translate our system into a vector-valued abstract Cauchy problem with suitable, non-standard 
coupled boundary conditions. The topic of partial differential equations on networks has become very 
popular in the last fifteen years, mostly due to its connections with quantum chaos which has motivated 
the introduction of the keyword "quantum graphs" . 

As natural for a theory tightly related to mathematical physics, partial differential equations associated 
with observables are usually considered in the theory of quantum graphs: so far, the large majority of 
investigations have been devoted to self-adjoint Laplacians and other second order elliptic operators. We 
mention the comprehensive survey |16j . On the other hand, also some ten years ago some pioneering 
investigations on system of transport equations on metric graphs have been commenced in |15) . The 
study of the first derivative on metric graphs inside the mathematical physical community dwells on the 
interpretation of i-4- as the momentum operator, and is much less common: we are only aware of [141 lllj . 
where some of the results of |15l [3] have been re-discovered and complemented by thorough spectral 
investigations. The aim of this note is to connect these two theories. The easiest case of the coupling of 
one diffusion and one transport equation (both possibly vector- valued) has been considered already in |12j . 
The work by Gastaldi and Quarteroni is presented by the authors as a first step towards the coupling of a 
Euler and a Navier-Stokes equation. A related viscosity analysis has been performed in jj. 

In |12j . well-posedness of the problem has been studied by means of viscosity methods under certain 
(relatively strict) coercivity assumptions. We are going to weaken the coercivity assumptions and study 
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a general system consisting of finitely many coupled intervals of either type. We are going to do so by 
considering a large class of coupled boundary conditions for operator matrices including both second and 
first derivatives. 

While the proposed dynamics actually seems to reflect the observed phenomena, it is very difficult to 
make an educated guess when it comes to propose natural transmission conditions. It seems that the search 
for correct transmission conditions, which is easy in the case of purely diffusive f |17|) or pure transport-like 
systems ( 15J), is much less trivial in the mixed case. As Gastaldi and Quarteroni put it in [12l § 1], "When 
coupling Euler and Navier-Stokes equations, the proper interface conditions are not obvious, in advance. 
A possible way of deducing them is to see the coupled problem as a limit of two coupled Navier-Stokes 
problems with vanishing viscous terms [in either one of the regions]". But even this limiting process is 
quite delicate, as the analysis in [T^l H] shows. 

In fact our main result, Theorem|4j states that systems of the kind discussed above are well-posed under 
a large class of transmission conditions. Thus, we complement a spectral analysis of the considered system 
by the study of typical features of semigroups for heat and transport equations, in particular dissipativity 
and positivity. This may hopefully suggest suitable conditions in specific settings, by allowing to reject 
large subclasses of boundary conditions as unfit. On the other hand, while many dissipative or dispersive 
extensions of A can actually be found, we are not sure that one can expect contractivity or positivity of 
the associated semigroup, when one thinks of some motivating systems. 

The plan of our article is as follows: In Section [2] we introduce the specific class of boundary conditions 
we are going to investigate. We discuss sufficient conditions in dependence of the boundary condition for 
well-posedness of the problem in the L 2 -setting and then in a general L p -setting, while in Section [3] our 
attention is devoted to the positivity properties of the solutions (observe that the space of L p -functions on 
a graph is a Banach lattice). In Section H we derive a secular equation for our system. 



Unlike in |12| , a possible motivation for the introduction of our setting arises from some biomathematical 

considerations. More precisely, it is known that electric signal coming from a neuron undergoes a certain 
synaptic delay before reaching another neuron, and cannot double back: this suggests to model this process 
by a system of diffusion (in the dendrites or axons) and transport (in the synapse) equations. In fact, it 
turns out that our setting can be adapted to discuss certain classes of coupled system of diffusion equations 
with boundary delays. This is explained in some detail in the conclusive Section [5] 

This article was partially completed during a visit of the first author at the University of Ulm under the 
financial support of the Land Baden- Wiirttemberg in the framework of the Juniorprofessorenprogramm - 
research project on "Symmetry methods in quantum graphs". 



We use the standard construction for the singular manifold on which the dynamic of the system is going 
to take place: We consider finitely many intervals of finite length that are connected to realize a simple 
finite metric graph G with node set V and edge set E. For the purposes of this note it is important to 
consider a partition of E in two disjoint subsets E^, E t , which are going to represent the edges on which 
diffusion and transport phenomena are going to take place, respectively. Hence, we denote by e^i , . . . , e^c 
and e tl , . . . ,e t ^ the edges of the metric graph. To each edge e^ or e t j we associate a length am or a t j, 
respectively, which in turn determines an orientation of each edge from to the other endpoint. More 
specifically, this orientation leads to the representation 



and in this case we write e(0) = v, e(a) = w. (Here and in the following a denotes the generic length - 
that is, either a<jj or a tJ -, depending on the context.) 

The structure of the network is given by the |V| x |E|- outgoing and ingoing incidence matrices T + := (1+) 
and T~ := (t^) defined by 




2. Boundary conditions for a mixed operator on a metric graph 



e = (v,w), 




for v, w e V, 



(1) 




and 
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The matrix I := (t ve ) defined by I := X + —T is the incidence matrix of G. Furthermore, let T(v) be the 
set of all edges incident in v, i.e., 

T(v) :={eeE: e(0) = v or e(l) = v} . 



For the sake of notational simplicity, if e = (v,w) we denote the value of a function u e : [0,a] — > C on e 
at and a by zt e (v) and w e (w), respectively. With an abuse of notation, we also set u e (v) = whenever 
e ^ T(v). By assumption, E is a hnite measure space. The space L 2 (E) of square integrable functions 
defined on the intervals associated with the edges in E becomes a Hilbert space with respect to the natural 
scalar product 



E 



T 
3=1 



Ue tj (s)v e t As)ds. 



By L 2 (Ed) and L 2 (E t ) we denote the L 2 (E)-functions with support on the edges belonging to E^ and 
E t , respectively. Finally, by H 1 (E) we denote the Hilbert space of square integrable functions that are 
weakly differentiable with weak derivatives in L 2 (E); and use the notation i/ 1 (E ( i) := H 1 (E) n L 2 (Ed) and 
H\E t ) :=H 1 (E)nL 2 (E t ). 

On L 2 (E) one considers the diagonal operator matrix 



(2) 



A := 



dx 2 







dx 



Au := A 


Ud 




'< 








-u' t 



We are going to present a class of m-dissipative realizations of A with domains containing the minimal 
one H^Ed) © Hg(E t ) and contained in the maximal one H 2 (Ed) © H 1 (E t ): that these are the extreme 
relevant cases follows from the observation that in both cases the transport and the diffusion part are 
clearly edge-wise decoupled. 

Remark 1. It has been observed in [16] that the most general class of self-adjoint extensions of the second 
derivative acting on the space of L 2 -functions over a metric graph with m edges ( each identified with an 
interval (0, 1) ) can be parametrized by a family of boundary conditions 



(3) 



P 1 



-u'(0) 
u'(l) 



(L + P) 



u(0) 
u(l) 



= 0, 



where P is any orthogonal projection of C 2m and L is any Hermitian matrix that can be seen as an 
operator on the range of P^ := 1 — P. (This kind of coupled conditions had already appeared in [7 in an 
investigation of Sturm- Liouville problems on graphs). Moreover, it has been implicitly observed in [9] that 
the first derivative acting again on L 2 (0, l) m 



(4) 



u(l) = Uu(0) 



generates a unitary group (i.e., by Stone's theorem, i-^ is a self-adjoint operator on L 2 (0, V) m ) if and only 
if U is a unitary m x m-matrix. When one tries to couple first and second derivatives by their boundary 
values, one soon observes that conditions of the form Q, while strongly motivated by the integration by 
parts formula, are not easily related to ^ - which in turn arise from the Gaufi-Green formulae. A possible 
way to overcome this problem is to write Q as 



(5) 



|u(l)| 2 = \Uu(0)\ 2 or rather Re(u(l) + Uu(0),u(l) - Uu(0)) = 0. 



We impose conditions that force the numerical range into the complex left half-plain, 
integration by parts yields 



An elementary 



(6) 



(Au, v) 



{u'd,v d ) - I {u t ,v t ) 

( U 'd: V 'd) + [u' d ,Vd]dE d 



+ / ( u t,Vf) - [u U V t ] d E t , 



E, 
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for u,v € H 2 (Ed) © H 1 (E t ), where here and in the following we denote 



[ud,v d ]dE d ■= ^ (udi(adi)vdi(adi) ~ u dl (0)v d i(0)) 



i=l 
T 



l u t,v t \ dEt 



^2 (uu(atj)vtj(a t j) - u tj (0)v tj (0)) 



3 = 1 

and for brevity |m( || Et := [ut,u t )o Et . While the equation ^ is not particularly appealing, the real part of 
the associated form is 

(7) Re(Au, u) = 



,' |2 



1 



Re[u' d ,u d ] aEd - -Mk> u G H2 ( £ d) x H\E t ). 



2' 

For the sake of notational simplicity we introduce the 2\Ed\ + |E t | dimensional ("boundary") Hilbert space 

U = U+ ®U d ®U t , -H*:=Cl Ed l, H t := C |Et| 
and define for all u g H 2 (E d ) x H 1 (E t ) the two vectors u,uGT-L, where 

{udi(a,di)}l<i<D 

(8) u:= {udi(0)}i<i<D and u 

2-i{utj(a tj ) + U tj (0)}i<j<T_ 



Wdi(adi)}i<t<D 

{-<i(o)h<i<D 

2~^{-u tj (a t j) + u t0 i(0)}i< j<t. 



are given with respect to the decomposition of W. With this notation, equation ^ becomes 



(9) 



Re (Au, u) 



.1 |2 



Re(u,u) n u e H 2 (E d ) x H 1 (E t ). 



These computations motivates us to introduce a class of boundary conditions of the form 

(10) P ± y L + (L + P)u = 0, 

where P is an orthogonal projector acting in H, P x := Id — P denotes the complementary orthogonal 
projector and the matrix L is an operator in the subspace KerPc % (whose extension by to the whole 
of % wc still denote by L). The boundary conditions can be equivalently written as 

P- 1 ^ + Lu = and Pu = 0. 

Finally we define the operator Ap^ which is studied in this work to be given by Ap^u :— Au on the 
domain 

D{Ap. L ) := {u€H 2 (E d )xH 1 (E t )\P x ^+(L + P)u = 0}. 

Remark 2. We call the boundary conditions ( |10[ ) type-decoupling if 

• P is an orthogonal projection of L 2 (E) onto {0}, W, C 2 ' Ed ' x {0}, or {0} x C' Et ' ; and additionally 

• L is a block-diagonal matrix with respect to the decomposition % = C 2 ' Ed ' © C' Et '. 

In other words, the boundary conditions are type- decoupling if actually no interaction between the boundary 
values in L 2 (Ed) and L 2 (E t ) takes place: this is of course the most trivial case, because the dynamics of 
the system can be effectively reduced to that of two distinct, non-interacting systems - a diffusive one and 
a transport one. 

Lemma 3. The vector space H 1 (E) := H l {Ed) ®H l (E t ) is a Hilbert space with respect to the natural inner 
product 

(u,v)hi(E) ■= {u\v') L 2 (£d) + (u',t/) L 2 (Et) + (u,v) L 2 {E} . 
The same holds for {u € -ff x (E) : Ud G H 2 (Ed)} = H 2 {Ed) © H 1 (E t ) with the inner product 

{u" ,v") L 2 {Ed) + (U, V) if 1(E)- 

Their embedding into L 2 (E) is dense and of p-th Schatten class for all p > 1. Furthermore, there exists 
C > such that the Gagliardo-Nirenberg-type estimate 

(11) \\u\\lo 0(E) <C\\u\\ L 2 iE) \\u'\\ LHE) for all u € H 1 (E) 
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holds. 

Proof. Throughout this note we are assuming our graph to be compact, hence the embedding of 
H 2 (E d ) © H l (E t ) into L 2 (E d ) © L 2 (E t ) is of p-th Schatten class for all p > 1, cf. [13]. 

A constant C > can always be found so that (11) holds, as this inequality can be reduced the case of 
intervals. Let namely u £ H 1 (Q, 1), then 

(12) Ki)l 2 <2V2|M| L 2 (0!l) |M| HW) . 

Indeed if u £ H 1 ^, 1) is such that tt(0) = 0, we can write 

\u{y)\ 2 < [ {u 2 )'{x)dx for all ye [0,1], 
Jo 

and by Cauchy-Schwarz's inequality we get 

\u{y)\ 2 < 2||u|| L 2 (0jl) ||ti / || i 2 (0)1) for all y e [0, 1]. 

For a general u £ i? 1 (0,a), it suffices to apply the previous estimate to u(x) = (1 — -)it(-) + -«(-). To 
derive ( 11 ) on the whole graph it suffices to sum up the estimates obtained on each edge separately. □ 

We are finally in the position to state our main result. 

Theorem 4. Let be P an orthogonal projector acting on % and L a linear operator on KerP. Then the 
operator Ap^ is quasi-m-dissipative. It is even m-dissipative whenever —L is dissipative. Accordingly, 
Ap,L generates a quasi-contractive (contractive, if —L is dissipative) semigroup on L 2 (E). 

Of course, the semigroup generated by Ap l is reducible (i.e., it leaves some non-trivial ideal of L 2 (E) 
invariant: in this case, some function space over a subset of E d or E t ) if the boundary conditions (10) are 
type-decoupling in the sense of Remark [2] 

Before proving this theorem, we need two preparatory lemmata. 

Lemma 5. Let u > be such that —L — uj is dissipative. Then 

u?C 2 



Re(Ap t LU,u) < — - — |MIl2(e) for u all satisfying (10) u£ D(Apj J ), 
where C > is a constant, depending only on the total length of the graph, such that 
(13) |u| 2 < C||u|| L 2 (E) ||u'|| L 2 (E) for allu£H 1 {E). 

Of course, such an u always exists: simply take uj := ||L||, although this is in general not optimal. 
Proof. Take u £ D{A PX ). Then it follows from ([9) that 

li 2 (E d ) 

If — L is dissipative, then we just estimate this by 



Re(A P , L u,u) = ||w'|||2 (E . ) + Re(Lu,u) n . 



-Re(Ap L u,u) > \\u'\\ 2 L2(E) > 0. 
If instead — L is only quasi-dissipative, i.e., if — L — uj is dissipative, then we deduce that for all e > 

- Re(A P , L u, u) > ||t/||2- 2(E) + Re(Lu,u) n 
(14) ^ l|w'|||2 (E) - wC||u'|| L 2 (E) ||u|| i 2 (E) 



In particular, for e — ^ we obtain 



> h'|l!2 (E) -^(f|| M '|| 2 L2(E) + i|| M ||| 2(E , 



^ 2 c 2 „ , |2 



-Re(A P . L u,u) > j-|M| L 2 (E) . 

This concludes the proof. □ 
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The restriction of A to H 2 (E d ) © #o(Et) is denoted by A . The adjoint of A in the Hilbert space L 2 (E) 
is the operator 



B := 



4 



£>(B) -P^eP, 1 ^ 



As Ap t L is an extension of B* , A* p L is a restriction of P. Therefore it can be described in terms of boundary 
conditions imposed on B. We introduce the notation 



{vdi{ai)}i<i<D 

{vdi(0)}l<i<D 

2~'{(vtj(0) +u*j(aj))}i<j<r. 



and v = 



{*40i)}l<i<-D 

{-'4(o)}i<<<d 

2- J {(-Wy(0)+Uy(a J -))}l< i <T. 



observing that 



f = v, and S = Ju, where J := 





-lci E * 



(The change of sign in the last component is due to the change of the direction on the transport edges). 

Lemma 6. The adjoint operator of Ap ^ is given the restriction of B onto 

D(A* PL ) :={(v d ,v t )eH 2 (E d )®H 1 (E t ) | (L* + P)v + P^ = 0}. 

Proof. By definition, the adjoint of Apx is the operator given by 

D{A* P L ) := {v G L 2 (E) | 3u G L 2 (E) s.t. (A P>L w, v) = (w, u) for all w G £>(Ap,z)}, 
A* PL v := u. 

To begin with, observe that following the computations in ^ the operator A - without boundary conditions 
- satisfies 



(Au,v) — (u,Bv) = — (u,v)<n + (u,v) 



H 



or rather 

(15) (Au, v) - (u, Bv) = 



Ik 




" 






v. 


U 




1« 







Jv 



u,v^H 2 {E d )xH 1 {E t ). 



Recall that A* p L is a restriction of the operator B. Hence v G P(Ap L ) if and only if the boundary term 
in (15) vanishes for all u G D(Apj J ). The range of 

[•]p.l : D(A P , L ) -> H 2 , [«] Pi£ = («,&) 

is Ker(L + P, P- 1 -) and hence the boundary term vanishes for a fixed u G H 2 (E d ) x P 1 (E t ) and all u G 
P(^4p.l) if and only if 

(«, Jv) JL Ker(-P- L , L + P), 

taking into account that 

° t J Ker(i + P, P- 1 ) = Kcr(-P\ L + P). 

Since the orthogonal complement of the space Ker(— P- 1 , L+P) is exactly Ker(L*+P, P- 1 ) one summarizes 
that v G L>(A PL ) if and only if <u G P 2 (E d ) x ff^Ej) and 

Lw) = (v, v) G Ker(L* +P, P^). 
This completes the proof. □ 
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Proof of Theorem^ It is known ([TU1 Cor. II. 3. 17]) that a sufficient condition for a densely defined operator 
to have m-dissipative closure is that both it and its adjoint are dissipative. By Lemma [5] Ap t L is quasi- 
dissipative for any L and dissipative whenever — L is dissipative. Like in Lemma[5]one proves that conversely 
the operator A* p L is quasi-dissipative for any choice of L and dissipative whenever — L is dissipative. 

To conclude the proof it suffices to check that Ap^ is actually closed. Because both the first and the 
second derivative without boundary conditions are closed operators, in our case it suffices to check that 
the boundary conditions are respected in the limit. This follows from the fact that the operators u n- u 
and u are bounded from H 2 (Ed) © H 1 (E t ) to %. □ 

Remark 7. Observe that dissipativity of the matrix —L is sufficient but not necessary for Ap^ to be 



Re(A P , L u,u) < -\\u'\\ 2 l2(e) +ojC -||u'||£a fE) + — \\u\\ 2 L2(E) for all u £ D{A PiL ) 



m-dissipative. By (14) one obtains the estimate 

L 2 (E) -tuj^ I ^IKIll^E) + ^11" L-iE 

for some C > only depending on the total length of the graph and all e > 0. 

To give a concrete example of an m-dissipative operator Ap^ where —L is not dissipative we consider 
the graph consisting of one transport and one diffusion edge with length a\u = ol\t = 1/2, for some I > 0. 
We consider boundary conditions as in ^ by taking 

"10 



P 



2- 1 2" 1 
2- 1 2" 1 



ami Lc ■- — —P^, 



c 



Re(A PtL u,u) < --||m / |||2 (E ) + ^||w|||2 (E) for all u G D{A P . Lc ). 



where C > is equal the constant appearing in the above estimate. We take advantage of the degrees of 
freedom that are left: Set ui := \\Lc\\ = C , and plugging e = 1 in the above inequality we obtain 

Since P defines a Dirichlet boundary condition on one endpoint of E^, by the one- dimensional Poincare 
inequality [31 Satz 5.16] there exists an a > 0, independent of I, such that^ 

2c? 

II M 'II^(E) > -jrNll»(E) f° r aU u e D ( A P,Lc)- 

Therefore if I > a 



Re(A P , L u,u) <l—--j \\u\\i HE) < 0, 
i.e., Ap^ c is dissipative even if the rank one operator ~Lc = qP 1 ' was not. 

Corollary 8. Let —L be dissipative. If additionally the matrix semigroup generated by —L on % is 
contractive with respect to the oo-norm, then the part of Ap ^ in L P (E) generates a semigroup of contractions 
in L°°(E) for each p £ [2,oo], which is strongly-continuous for p < oo. 

If the matrix semigroup generated by L* on % is contractive with respect to the oo-norm, then the 
semigroup generated by Ap^ in L 2 (E) extrapolates to a strongly- continuous semigroup of contractions in 
LP{E) for eachp£ [1,2). 

Recall that the semigroup (e tA ) t >o generated by an n x n matrix K — (k^) with complex-valued 
coefficients is oo-contractive if and only if 

(16) Reku + \kih\ < for all i = 1, . . . , n, 



1 Observe that functions in D(Apj JC ) are in general not of class H 2 (0,l) (in fact, not even H 1 (0,i)), but this is not 
relevant: in order to apply Poincare's inequality it suffices to have a function that is a.e. weakly differentiable, and the 
possibilty to find a path from each point of the graph to the endpoint(s) where the function vanishes; and that the boundary 
values are uniquely determined along this path - this follows in our case from the fact that ut(0) and ut(at) are uniquely 
determined, since the system 

ut (0) + u t (a t ) = V2u d (0) and - u t (0) + u t (a t ) = V2u' d (0) 
admits a unique solution for all u d (0),u' d (0). 
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cf. [TS1 Lemma 6.1]. 

The proof of this result is based on a generalization due to Yokota of Brezis' classical results for invariance 
of closed convex sets of a Hilbert space under the semigroup generated by a subdifferential. In the linear 
case, his results read as follows, cf. [351 Thm. 2.4]. 

Lemma 9. Let X be a complex Hilbert space and K C X be closed and convex, and denote by Pk the 
orthogonal projector onto K . Let T — u) be an m-dissipative operator for some oj > and (S(t)) t >o be the 
strongly-continuous semigroup generated by T . Then K is invariant under (S(t)) t >o if and only if 

Re(Tu, u - P K u) < u\\u - P K u\\ 2 for all u £ D(T). 

Proof of Corollary^ Let K be the unit ball of L°°(E). The associated orthogonal projector is given by 

Pk u '■— (\u\ A 1) sgnM and therefore u — Pku :— (\u\ — 1) + sgntt, «eL 2 (E). 

Accordingly, 

(P K u)' := u'l{ u <i} and (u - P K u)' := u'l {u > 1} , ueH 1 (E). 

Let us study the condition in Lemma [9j which in our case reduces to 

Re(A P ^ L u, u - P K u) < for all u £ D(A PX ). 

Writing with an abuse of notation Pku = (PkUcI, PRUt) (this is not completely inconsistent, due to locality 
of Pk), one has for all u € H 1 (E) 



Re 



{u'd,(u d - PKU d )) = - Rc[u d , {u d - PKU d )')dE d + Re / {u' d , (u d - PKUd)') 
= - Re[u' d , (u d - PKU d )']dE d 

+ Re / ((P K u d y, (u d - P K Ud)') + Re / {{u d - P K u d )' , {u d - P K u d )') 



Re[u' d , {u d - PKU d )'}dE d + I \\{ud ~ PKUd)'\\ 2 

and likewise 

Re/ (u' t , (u t - P K u t )) = Re I ((P K u t )', (u t - P K u t )) + Re / ((u t - P K u t )' , (u t ~ P K u t )) 



whence 



Re/ \\{u t -P K u t )\\ 2 -Re [ (u' t , {u t - P K u t )) , 



Re jf (u' t , (u t - P K u t )) = \f^ II K - PKUt)\\ 



(Here each disappearance of an integral terms follows from the fact that the integrand is the product of 
two functions with disjoint support). We conclude that for all u € D(Ap^) 

Re(A P . L u,u- P K u) = - [ \\ (u d - P K Ud)'\\ 2 - Re(L(u - P K u), (u - P K u))u < 0, 

J Ed 

where the last inequality follows from dissipativity of —L. Hence, under our assumptions the semigroup 
generated by Apx is both L 2 (E)- and L°°(E)-contractive. By the Riesz-Thorin interpolation theorem, 
each semigroup operator e tAp - L extends to a bounded (in fact, contractive) linear operator on L P (E), for 
all p £ [2,oo]. Strong continuity for all p < oo is also clear. The second assertion can be proved by 
duality, taking into account Lemma [6j Strong continuity of the semigroup in L 1 (E) is not obvious, but 
is a consequence of [131 Prop. 5]. In fact, by Lemma 12 below Ap : L has compact resolvent and therefore 
by [IHl Cor. 2.15] each orbit of the generated semigroup is relatively compact, hence weakly relatively 
compact. □ 
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Remark 10. Instead of considering the transport- diffusion-type Cauchy problem associated with dissipative 
operators, one might think of a Schrddinger-type Cauchy problem involving self-adjoint operators with mixed 
dynamics. To this aim, consider the symmetric operator 



5° := 



d?_ 

dx 2 





Q 

dx 



D(S°):=H^(E d )®H^E t ). 

It has equal deficiency indices (2| j + \E t \, 2\E d \ + |E t |) and its adjoint S — (S )* is formally the same oper- 
ator with domain H 2 {E d )®H l (E t ). Hence there exists self-adjoint extensions and these can be parametrized 
in terms of boundary conditions. One defines the appropriately modified vectors of boundary values 



{ u dii.O>di)}l<i<D 
{udi(0)}l<i<D 

2 {u t j(a tJ ) +Utj(0)}i<j<T. 



Wdi( a di)}l<i<D 

{-<i(Q)h<i<n 
~ 2 {-u t] {a t] ) +u tJ (0)}i<j<T. 



to obtain the well known Hermite symplectic form on the space of boundary values 



(Su, v) — (u, Sv) 



u 




" 






V 


u 











V 



u,v £ H 2 (E d ) © iJ 1 (E t ). 



'hi- 



lt is known that there is a one-to-one correspondence between the self-adjoint extensions of symmetric 
operators and the maximal isotropic subspaces with respect to this Hermite symplectic form. A unique 
parametrization of such a subspace is given in terms of a projection P and a Hermitian operator L acting 
in KerP . Therefore any self-adjoint realization Sp : L of S° is a restriction of S onto a domain of the type 

D(S P , L ) = {u£ H 2 (E d ) ® H\E t ) \Pu = and Lu + P^u = 0}. 

The spectrum of the self-adjoint operator Sp^ is purely discrete and there are sequences of eigenvalues 
going to +oo as well as to — oo. These self-adjoint operators can be interpreted as Hamiltonian consisting 
of a standard Laplacian and a less usual moment-type observable. As already mentioned, moment operators 
on graphs have been recently studied in [11] . By Stone's theorem, for self-adjoint Sp^L the Cauchy problem 



'dt 



U(;0) 



is governed by a unitary group. 



Sp,l) u(x,t) = 0, 
= u 

3. POSITIVITY 



t £ 



Observe that L 2 (E) is a Hilbert lattice, hence we can discuss positivity of the semigroup. 

Theorem 11. Let P an orthogonal projector in % and L a (2 1 EE,^ | x |E t |) x ( 2 1 EE^ | x \E t \)-matrix such that 
L = P ± LP' L , where P = Id — P. Then the following assertions hold. 

(a) The semigroup generated by Ap^ leaves the set of real-valued functions invariant if and only if L has 
real entries. 

(b) The semigroup generated by Ap ^ is positive if and only if L has real entries that are positive off- 
diagonal. 

Proof, (a) By Lemma [9] and reasoning as in the proof of [201 Prop. 2.5], one sees that the semigroup 
generated by Ap^ is real if and only if {Ap^ Re u, Imu) £ R for all u £ D{Apj J ). Taking into account (J6|, 
a direct computation shows the assertion. 

(b) By Proposition [5j Ap^ is always quasi-dissipative, say with constant lu. By [THl Thm. CTI.1.2] it 
suffices to show that our assumptions on L are satisfied if and only if Apx — w is dispersive. However, 
because the semigroup generated by Ap^ — uj is clearly positive if and only if so is the semigroup generated 
by Ap^, we can confine ourselves to the case of Ap^ dissipative. In other words, it suffices to check that 
our assumptions on L are satisfied if and only if for every u £ D(Apj j ) - in fact, by (a), if for every 
real-valued u £ D{Apif) - there exists < (f> u £ L 2 (E) such that 

• ||^u|U=(E) < 1, 

• (u, <Pu)l 2 (E) = \\u + \\ L i(E) and 

• (Ap LU, 4>u)l 2 (e) < 0. 
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Let in fact u € D(Ap^) be real- valued and define <f> u as a vector in the positive cone of L 2 (E) by 



Then, one has 



Furthermore, 



^lll 2 (E) 



1 



li=(E) 
f 



u+| 


L 2 (E) 


u 


t 


u+\ 


L 2 (E)_ 






L 









,+ |2 



= 1. 



II u MIl 2 (e) 
1 

l|w + IU 2 (E) 
l w+ IU 2 (E) 



,+ |2 



Et 



,+ |2 



E f 



Finally, we take a closer look at ([6]) and observe that 

(u' t ,u t ) l{ Ut >0} = 



E t 



E t 



-<n t ,^}+ / (u t ,(u+)') 

-(u t ,uf)dE t + (Ut,u' t )t{u t >0}, 



E, 



whence 



i + \\l 2 (e)(Ap,lu, 4> u 



(ulu+)- / (u' t ,v t ) 



E t 



ll u dllL(E d ) 1 {« £i >0} + (w'd^rf)^ - -(u u ut)g Et 
\\ u 'd\\L2(E d ) 1 {u d >0} + (^,W + )c 2 l E £ ilxl E t 
ll M 'dlli2(E d ) 1 {«d>0} - (iM,M + ) c 2|E d |x|E t |. 

Since the boundary values of u and u' that appear in u are mutually independent, one concludes that 
{Ap^u, 4> u ) < if and only if 

(17) (i$,e + > R2 |E d ix,E tl > foralUe R 2|Ed|x|Etl . 

Because £ + is the orthogonal projection of £ onto the positive cone of R 2 l Ed l x l E *l and (L-, •) is the quadratic 



form associated with L, we deduce from the invariance criterion [201 Thm. 2.2] that (17) is satisfied if and 
only if the semigroup generated by L is positive, i.e., if and only if L has real entries that are positive 
off-diagonal. □ 

We close this section by formulating a conjecture on the general behavior of the semigroups generated 
by the m-dissipative operators Ap^. 

Conjecture 1. Let E f 7^ 0. If a function f is supported in e € E t , then the semigroup generated by Ap ^ 
will shift its profile until the support of e tAp L f hits an endpoint of e. Because e tAp - L f has in this lapse 
of time the same profile of f, the semigroup cannot be immediately smoothing. In particular, it cannot be 
analytic - in fact, not even immediately differentiable. On the other hand, it seems reasonable to imagine 
that the semigroup smoothens the profile of a function as soon as it reaches an edge in E^. We conjecture 
that if Ed 7^ 0, then the semigroup is differentiable for all t > T , where T is the length of the longest path 
inside E t (to compute taking into account the possibly coupled boundary conditions). 
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4. Spectral theory 

Taking into account Lemma [3j we promptly obtain the following. 

Lemma 12. For all orthogonal projectors P inJi and all linear operators L acting in KerP, the operators 
Ap,L have resolvent of p-th Schatten class for all p > 1 (and, in particular, of Hilbert- Schmidt class). In 
particular, Ap^ has only pure point spectrum. 

4.1. Non-zero eigenvalues. In order to determine the pure point spectrum of Apx, a natural Ansatz 
for finding eigenfunctions is to take k E C \ {0} and to consider 



a di {h)e ik * + f3 di (k)e~ ikx , 
ltj(k)e k2x , 



</>{x,k) = 

The boundary conditions (P + L)<p(-, k) + P J -0(-, k) = are encoded into 



X (E ^di 7 ^ L 

x e e t j, j = 1, 



,D, 



where {a d (k)}i-i t ... lD = 
coefficients. The matrices 



[{P + L)X{k)+P^Y{k)\ 
a d i(k), {/9d(fe)}i=i,...,D = Pdi{k), {lt(k)}i =] , 



a d (k) 
Mk) 
7t(k) 



,t = Jtj(k) are the sought after 



X(k) 



Y(k) = 



1 



—ik 




V2 

ike~ ik - 
ik 







Wi + e fc2 a t ) 



act in H are given with respect to the decomposition H = T-CJ © W d © Ht ■ Here the notation 



!,.; , / = 5 ie ±k ' at - and {< , D = 6 3l e ±ika * 



is used. Accordingly, the following holds. 

Proposition 13. For all orthogonal projectors P in Ji and all linear operators L acting in KerP , the 
number —k 2 € C \ {0} is an eigenvalue of Ap r, if and only if the matrix 

Z(k) := [(P + L)X(k) + P X Y(k)}, 

has non trivial null space. The geometric multiplicity of —k 2 equals the dimension of KerZ(k). 

Hence, the secular equation is 

detZ(fc) = 0. 

Example 14. Consider the graph consisting of one diffusion edge of length ao and one transport edge of 
length ap- Let P be as in Remark^ 

"10 " 

2- 1 2- 1 
2- 1 2- 1 



P 



and L = 0. 



Then the operator Ap l is m-dissipative and the secular equation becomes 



det Z(k) 



V2 



an(fca£>) ^1 



k COS 



i(ka D ) ( 



1 + e 



k a t 



0. 



In particular the spectrum of Ap^ contains a sequence of real eigenvalues going to — oo. 

It seems in general difficult to give precise statements on the distribution of the eigenvalues. 
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4.2. Eigenvalue zero. For the eigenvalue zero one uses for the eigenfunctions the Ansatz 



(j>o{x) 



otdi(0) + Pdi(0)x, x e e dl , i = l, 
7tj(0), x e e tf , j = 1, 



,D, 
.,T. 



The boundary conditions (P + L)<f>(-, 0) + P^4>(-, 0) = are encoded into 



[(£ + P)X + P^Fq] 



'ad(O)' 
&(0) 
7t(0) 



with 



1 a d ' 
10 
V2 



and Yr 



o 



10 
0-10 




where {&,j}j,i=x t ...,D = Sjia-dj- Hence we obtain a characterization for the eigenvalue zero. 

Proposition 15. For all orthogonal projectors P in H and all linear operators L acting in KerP, the 
operator Ap^ is invertible if and only if 

detZ o ^0, where Z Q = (P + L)X Q + P^Yq. 

In particular the invertibility of the operator Ap t L is independent of the lengths of the transport edges. 

Example 16. Consider again the graph given in Example\T^ and let P be the projector given in Example \lJ\ 
and Remark^^ Let be L(j — CP 1 - , where P 1 - — Id — P and for arbitrary C G C Then 

detZ = -2"s - Ca D 23 

and therefore Apj J(J is invertible only for C ^ — (2au) . 

4.3. The resolvent operator. Knowing the resolvent R(X, Ap^) = (Ap.L — A) -1 for Ap^ dissipative 
and for A in the right halfplane, one can describe the semigroup generated by Ap ^ by means of the inverse 
Laplace transform as 



(18) 



e tAp - L u 



= lim / e tx R(\,A PL )ud\, ueL 2 {E d )xL< 
J e _ in 



(Et), 



for any e > 0, see [H Thm. 3.12.2]. In fact, in the following an explicit formula for the resolvent is given 
in terms of the boundary conditions and the edge lengths. First we define the shorthand notation 



E 



/ u t j(x)dx, for u € L 2 {Q). 



(i 



Proposition 17. For all orthogonal projectors P in %, all linear operators L acting in KerP and for all 
k 7^ such that —k 2 € p(Ap^if), the resolvent operator R(k) — (^4p,L + k 2 )" 1 is the integral operator given 
by 

R(k)u(x) :— J r(x,-,k)u(-) 

with kernel 

r(x, y, k) := {r (x, y, k) - $(x, fc)£(P, L, fe)*(y, k)} W(k). 

Here we have denoted 

S(P, L, k) := [(P + L)X(k) + P ± F(fc)]- 1 [P ± P 1 (fc) + (L + P)R 2 (k)} 
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and furthermore 



r (x,y,k) 
W{k) 

{rd(x,y,k)} n ,m 

{rt(x,y,k)} jt i 





rd{x,y,k) 







r t (x 












lci E *i_ 




— X Jk\x-y\ 


m = 1, 


/ ( k 2 {x-y) 

= Sn r 


a; < y 




VI ' 


a; > y / 



, i,/ = l,...,T. 



Finally, 



$(ar, k) := 















$>(x,k) := 



Jky 
— iky 








-fc 2 y 



w/iere i/ie entries are diagonal matrices whose entries are functions with arguments from the corresponding 
edges and 



R^k) :-- 



'ikl c \E d \ 







ike tk - d 









ifc(fc) := 













Proof. It is sufficient to prove that r(x, y, k) is the Green's function of the operator (Ap^ + k 2 ). Consider 
the unperturbed operator 



R (k)u(x) = r (x, ;k)W(k)u(-) for u e C o °°([0, a*]) C o °°([0, ay])J . 

The equation (A + k 2 )R Q (k)u = u is satisfied on the diffusion edges as ^ e lk \ x -y\ i s the Green's function of 
Ld(k) = ^2 — k 2 on the whole real line (this follows from standard arguments using the Fourier transform). 
By continuing functions Uj G L 2 ([0, ct^i]) trivially onto the real line the claim follows. Similarly the diagonal 
entries of r t (x,y,k) are the Green's function for L t (k) = — ^ — k 2 on the whole real line, which follows 
from standard arguments from the theory of ordinary differential equations. Again by continuing functions 
Uj in L 2 ([0 7 a t j}) trivially onto the real line the claim follows. For the correction term one has 



(A + k 2 ) J $(x,k)E(P,L,k)V{;k)u{-)=0. 



Therefore (A + k 2 )R{k)u = u. As (0^ C§°([0, a di ]) 0^ =1 C§°([0, ay]) J is dense in L 2 (G) and r(-,-,k) 

defines for all k ^ such that ~k 2 e p(Ap^) a bounded linear operator on L 2 (G x G). One concludes by 
density that (A + k 2 )R{k)u = u for all u e £ 2 (G). 

It remains to prove that R{k)u € D(Ap,i)- Observe that for all a > and all u e L 2 (0, a) 



ifc|x— j/| 



u(y)dy 



d_ 

dx 
d 
dx 



f e lk \ x -y\u{y)dy 
Jo 

f e tk \ x - v \u{y)dy 
Jo 

f e^ x -^u{y)dy 
Jo 



c=0 JO 

=e 



x=0 



e lky u(y)d yi 

ka P e~ ik «u(y)dy, 
Jo 

ik \ e lky u(y)dy, 
Jo 



-ike 



ka I 
JO 



-iky 



u{y)dy, 
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and considering only the edge [0, a] 



r t (x,y,k 2 )u(y)dy 



r t (x 7 y,k 2 )u(y)dy 



x=0 



e- k v f(y)dy, 



This gives in the matrix notation for u £ L 2 (G) and v G TL 
Ro(k)u(x) = R 2 {k) ( V{k,-)W(k)u(-), 



Therefore 



Ra(k)u(x) = Ri(k) j V(k,-)W(k)u(-), 



$(x, k)v = X(k)v, 
$(x, k)v = Y(k)v. 



r(x,;k)u(-)={R 2 (k)-X{k)X(P,L,k)) / *(-,k)W{k)u(-), 
g Jg 



r(x,;k)u(-)=(R 1 (k)-Y(k)E(P,L > k)) / V(;k)W(k)u(-), 

G J G 



and hence for all u € L 2 (G) 



(P + L) r(x, ; k)f(-) + P^ / r(x, -,*)/(■) = 0. 

JG JG 



This proves that r(-, •, k) is the Green's function for the mixed problem. □ 

Remark 18. The statement in Theorem^ (a) can also be proved also by observing that the kernel ?"(•,•, in), 
k > has real coefficients whenever the operator L have real entries. Therefore the resolvents (Ap t L~ n 2 )~ l 
map real-valued functions to real-valued functions for any k 2 > 0. Applying the inverse Laplace transform 



described in (18), one concludes that the semigroup generated by Ap^ is real whenever —L is dissipative 
and real. 



5. Delayed diffusion equations via mixed PDEs on graphs 

We conclude this note discussing a possible application of our theory to delayed diffusion equations. It 
is a classical idea, thoroughly developed e.g. in [5], that delays can be mathematically modeled by means 
of (auxiliary) transport phenomena. While this theory is classical whenever the delay is "distributed" - 
i.e., it acts on each point of the domain of diffusion -, it is less standard and technically more involved 
whenever the delays only concern the boundary values. The situation we want to describe is the following: 
We consider a ramified structure with an ongoing diffusion process. Incoming particles are absorbed in 
some nodes and stored there for some time (which depends solely on the features of the node itself), before 
triggering a flow in the incident edges. This phenomenon can be described by attaching to each such 
node a fictitious loop, which the particles have to cross before reaching the adjacent edges. A comparable 
approach has been presented in the recent article [5]. 

As a tentative motivation of the investigation of this class of problems, we briefly discuss a (much sim- 
plified) model of a dendrodendritic chemical synapse, referring the reader to any introductory textbook on 
neuronal modeling (e.g., [351 Chapters 2 and 9] and jH Chapter 5]) for basic notions and some miscellaneous 
models, including delayed ones. Also for the sake of simplicity, we drop the absorption term and therefore 
replace the cable equation usually discussed in theoretical neuroscience by a simpler diffusion equation. 
(This in fact a bounded perturbation, hence neglectable when it comes to discussing well-posedness) . We 
also remark that our model essentially applies to the case of less exotic axodendritic synapses, if we linearize 
the active transport phenomena typical of the latter. 

Consider two dendrites (modeled by two intervals ei,e2) that are incident in the synapse v, which is 
terminal endpoint of ei and initial endpoint of &2- The synaptic input coming from ei undergoes a delay 
of Tdci = 1 before reaching e 2 and cannot double back. For the sake of simplicity, we also impose sealed 
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end conditions on the first dendrite ei as well as on the second endpoint of e2, but of course longer chains 
of neurons may be modeled in a similar way. 

The synaptic input is an action potential that lets neurotransmitters be released by synaptic vesicles. 
Experimental observations suggest that no obvious (linear) algebraic relation exists between the pre- and 
post-synaptic potential in the dendrites - i.e., between the boundary values of the unknowns Ui,u 2 in the 
diffusion equations. We propose to discuss this system by a network diffusion problem with boundary delay 



(BD) 



ui{t,x) 


= u'{{t,x), 


t > 0, x £ (0,1), 


u 2 (t,x) 


= U 2 {t,x), 


t > 0, X £ (0,1), 


Ui(t,l) 




t > 0, 


-u' 2 (t,0) 




t > 0, 


«'i(*,0) 


= o, 


t > 0, 


u' 2 {t,l) 


= o, 


t > 0, 


ux(0,x) 


= fi(0,x), 


X £ (0,1), 


u 2 (0,x) 


= f 2 (0,x), 


X£ (0,1), 


Ui(t,l) 


= fde\(t), 


t£ [-1,0]. 



The first boundary condition, of Robin-type, can be interpreted by saying that part of the reaching the 
presynaptic nerve terminal is reflected into the entrance dendrite. The second condition is actually the 
relevant one: It reflects the fact that the conduction speed of the signal in the dendrites is approximatively 
constant, at least within the same cortical area. In other words, even if the postsynaptic potential will in 
general have a different amplitude from the presynaptic one, their speed of propagation will be the same. 
Finally, it is easy to convince ourselves that the above problem is undetermined if the last initial condition 
on the delay term is not imposed. 

In order to implement the delay feature into the node conditions, we introduce an edge edei and an 
unknown Udci to model the transport of neurotransmitters in the synaptic cleft between the pre- and the 
postsynaptic neurons, hence effectively introducing a delay phenomenon. 

Our aim is now to replace the fourth, delayed equation in (BD) by two node conditions in the endpoints 
of edei- More precisely, we impose that 



u[(t, 1) = u de i(i,0), 



t > 0, 



as well as 

-l4(t,0) =Udel(t,l), *>0. 

In other words, the flow of postsynaptic potential and the flow of presynaptic potential agree, even if its 
transmission undergoes a certain delay (which we have normalized). 

We are hence led to consider an (undelayed) initial boundary value problem 



(BD') 



ui(t,x) 


= u'{(t,x), 


t>0, X £ 


(0,1), 


u dcl (t,x) 


= ""del (*>&). 


t > 0, x £ 


(0,1), 


u 2 {t,x) 


= u 2 {t,x), 


t > 0, x £ 


(0,1), 




= U de l(i,0), 


t > 0, 




= 


t > 0, 




-u' 2 (t,0) 


= U de \(t, 1), 


t > 0, 




«i(*.o) 

u' 2 (t,l) 


= o, 


t > 0, 




= o, 


t > 0, 




ui(0,x) 


= /i(0,x), 


are (0,1), 




u 2 {0,x) 


= h(0,x), 


are (0,1), 




k u dc i(0,x) 


= /dd(0,x), 


are (0,1), 





i.e., we have got rid of the boundary delay by passing to the a larger state space. Observe that the above 
model is intrinsically non-symmetric in the sense that potential can only flow from dendrite ei to e 2 , but 
not vice versa. This is a typical feature of chemical synapses, as opposed to electric ones. 

One can check that the problems (BD) and (BD') are equivalent. It should be observed that transport- 
based synaptic transmission models are not very common in the literature. Due to their numerical and 
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analytic complexity they are actually often replaced by convective-diffusive (or even purely diffusive) mod- 
els. A convincing pleading of a transport approach can be found in |24[ § 2 and § 6]. 

To discuss the problem (BD') in our setting, observe that (BD') can be seen as an abstract Cauchy 
problem on L 2 (E) with boundary conditions as in (10), where jE^j =2 and |E t | = 1 with 



L := 



1 












































1 











-V2 


v/2 











1 



and P := 0. 



A direct application of Theorems [4] and 11 yields finally the following 



Proposition 19. The initial-boundary value problem (BD') is governed by a strongly- continuous semigroup 
(T(i)) t > on L 2 (E) — L 2 (0, 1) x L 2 (0, 1) x L 2 (0, 1). This semigroup is real but non-positive, i.e., there exists 
to such that T(to) is not a positive operator. 

Remark 20. Observe that the above result does not really depend on our choice to consider a delay interval 
edei of unit length: we may in fact replace (0, 1) by an interval of arbitrary length. Furthermore, explicite 
computations show that 

det Z = 

and therefore Ap l is not invertible for any edge lengths. Moreover, —L is not dissipative, hence Theorem^ 
cannot be applied to deduce contractivity of the semigroup that governs the problem. 

An explanation for the non-positivity of (T(t))t>o comes directly from the biological description of the 
considered model: Due to lack of physical contact between the presynaptical dendrite and the postsynap- 
tical one, the potential is typically not conserved. This phenomenon has already been observed in [121 § 2] 
for different transmission conditions. 



Remark 21. More generaly one can consider couplings between the transport edge and the two diffusion 
edges of the form 

u[(t,l) = u d ci(i,0), <>0, 
u[(t,l) = 7Ui(i, 1), t > 0, 

-u' 2 (t,0) = Su 2 (t,0)-u de i(t,l), t>0. 
For any choice ofj, S € R this defines a well-posed problem, since Theorem^ applies to the system. 

We conclude by observing that another possibility - perhaps more natural in certain different contexts 
- could be to impose transmission of either boundary values at both endpoints. It seems that this setting 
is not covered by our setting, but we are not able to guess whether this is a hint that such a system is not 
well-posed. 
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